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Abstract 

In this paper we investigate the short-time decoherence of single Josephson charge qubit (JCQ). The measure of decoherence 
is chosen as the maximum norm of the deviation density operator. It is shown that when the temperature low enough (for 
example T — 30mK), within the elementary gate-operation time r 9 ~ 12.7ps, the decoherence is smaller than 10 -4 at present 
setup of JCQ. The Josephson charge qubit is suitable to take the blocks for quantum computations according to the DiVincenzo 
low decoherence criterion. 
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I. INTRODUCTION 



David DiVincenzo put forward five criteria for the can- 
didates of quantum computing hardware to be satisfied 
0. One of which is the low decoherence. An approxi- 
mate benchmark of the criterion is a fidelity loss no more 
than ~ 10~ 4 per elementary quantum gate operation. 
The Josephson junction is considered to be a promising 
physical realization of qubit. So it is a very significa- 
tive work to investigate the decoherence of the qubits 
based on the Josephson junctions. To perform quanti- 
tative studying of the decoherence for a given system, 
in general, one must firstly solve the quantum dynamics 
problem of the system coupled to the environment. Many 
kinds of methods have been used for this purposes Q 
The decoherence of the qubit based on Josephson junc- 
tions have been studied by many authors. In their inves- 
tigations the dynamics appeal to the Redfield formalism 
[4J or Bloch-type [5| master equations, where the Markov 
approximation are used [(| . The approximation is usually 
used to evaluate approach to the thermal state at large 
times. However, a quantum gate operation is finished in 
a instantaneous time so the approximation may be not 
suitable to investigate the decoherence of the qubit gates 
for quantum computing purposes. Recently, V. Privman 
et al. introduced two measures to quantify the short- 
time decoherence Q [3 ■ The measures are based on the 
short-time approximation of the split-operator and de- 
rived from a operator norm ||^4|| (seeing the following 
definition). By using these measures V. Privman et al. 
investigated some spin-boson models. In this paper we 
use one of the measures investigating the short-time de- 
coherence of the single Josephson charge qubit (JCQ) 
operations. 



paper. The single JCQ Hamiltonian is Q 



Ech (n 



Ej cos ip. 



(1) 



Here, E ch = e 2 / (C g + Cj) is the charging energy; Ej = 
I c H/2e is the Josephson coupling energy ,9], where I c 
is the critical current of the Josephson junction, h the 
Planck's constant divided 2ir, and e the charge of elec- 
tron; n g — C g Vg/2e is the dimensionless gate charge, 
where C g is the gate capacitance, V g the controllable gate 
voltage. The number operator n of (excess) Cooper pair 
on the island, and the phase ip of the superconducting or- 
der parameter, are quantum mechanically conjugate. Be- 
cause the Josephson coupling energy Ej is much smaller 
than the charging energy E c h, and both of them are much 
smaller than the superconducting energy gap A, the 
Hamiltonian Eq. Q can be parameterized by the num- 
ber of the Cooper pairs n on the island as j5| [Hj 



Hr = \ E ch (n - n g ) 2 \n) (n\ 



~Ej[\n) ( n +l\ + \n+l) (n\ 



(2) 



When n g is modulated to a half-integer, say n g = 1/2 
the levels of two adjacent states are close to each other, 
the Josephson tunneling mixes them strongly. When the 
temperature T is low enough the system can be reduced 
to a two-state system (qubit) because all other charge 
states have much higher energy and can be neglected. So 
the Hamiltonian of the system can approximately reads 



H s — — -B z a z — -B x a Xl 



(3) 



II. JCQ HAMILTONIAN AND A MEASURE OF 
DECOHERENCE 



In this section we firstly review the JCQ model and 
then introduce the measure of decoherence used in this 



where B z — E c h (1 — 2n g ) and B x = Ej. Eq.© is sim- 
ilar to the ideal single qubit model |5j, but it can be 
modulated only by changing one parameter B z . How- 
ever, changing the parameter B z (through switching the 
gate voltage) one can perform the required one-bit op- 
erations. If, for example, one chooses the idle state far 
to the left from the degeneracy point, the eigenstate lose 
to |0) and |1) . Then switching the system suddenly to 



1 



the degeneracy point for a time r and back produces a 
rotation in spin space p|, 



Uj — exp 



iEjT 



cos^- isin 2 ^ 
isin4^ cos^fi 



(4) 



This can be obtained by modulating the gate voltage and 
making B z — E c h (1 — 2n g ) = 0. So it is interesting for 
our to investigate the decoherence of this kind of qubit. 
If we consider the interaction of the qubit and its envi- 
ronment the total Hamiltonian becomes 



H = H X + H,+ H f 



(5) 



where Hb is the environment Hamiltonian which is usu- 
ally modelled by a bath of an infinite number of the 
harmonic oscillator models which is equivalence to an 
infinity-mode electromagnetic field. The dissipation and 
the decoherence of the quantum systems is considered 
because of the coupling of the quantum system and the 
fluctuating electromagnetic field. In the JCQ model, the 
coupling of the qubit with the electromagnetic fluctua- 
tions can be modeled by a impedance Z (to) , placed in 
series with the voltage (see Fig.l of Ref.[l2j). The bath 
Hamiltonian Hb and the interaction Hamiltonian Hj are 



Hi = A. 



(6) 



where A s = a z and 

M fe = Lu k a\a k , N k = g k a\ + g* k a k . (7) 

Here, oj k are the bath mode frequencies, a k and a k are 
bosonic create and annihilation operators. N k are the 
freedom of environment and it is direct proportion to the 
fluctuations of the voltage of the external circuits [f| \H\ 
and g k are the coupling constants. 

Because the bath modes are infinite, their frequencies 
can be taken continuous. The spectral density of the 
continuous bath modes is 



J v (w) 



k 



g\b (w - uj k ) 



(8) 



because it can result in a same Johnson-Nyquist relation 
which can be obtained from fluctuation-dissipation theo- 
rem on this system. On the other hand, at low-frequency, 
the spectral density behavior is 



J v (lo) — rjoj s exp (lo/lu c ) . 



(9) 



where uj c is the high-frequency cut-off and rj is the dimen- 
sionless strength of the dissipation. According to Ref . [T^| 
we know that in the particular interest ohmic case (s = 1) 
[l3| , for the model of JCQ coupling with electromagnetic 
fluctuations (see Fig.l of Ref. 

- A 

^ Ro 
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E2), 



Ct_ 
Cj 



(10) 



where Rq = (2e) /h is the (superconducting) resistance 
quantum. From Eqs.JSJ) and J§J| we know that for a ohmic 
bath 



glW{uj) = yywexp (lu/uj c ) , 



(11) 



where W (u>) is the density of states. 

There are many measures to characterize the decoher- 
ence. Usually the environment being assumed to be a 
large macroscopic the interaction with it leads to the 
thermal equilibrium at temperature T. In this case, 
Markovian type approximations can be used to quan- 
tified the decoherent process and it usually yields the 
exponential decay of the density matrix elements in the 
energy basis of the Hamiltonian H s . In this time scale 
the measures of entropy and the first entropy are used 
for quantifying the decoherence. But the decoherence 
of the qubit gate operations cannot be characterized by 
this methods because the time of the elementary quan- 
tum gate operation is much shorter than the thermal 
relaxation time. It has been shown that the norms ||c|U 
is useful for describing the decoherence of the short-time 
evolution. Here a is the deviation operator defined as 



f (r) = p (r) - p* (r), 



(121 



where p (r) and p 1 (r) are density matrixes of the "real" 
evolution (with interaction) and the "ideal" one (without 
interaction) of the investigated system. |jcr|j A is defined 
as 



( i<p\<r\<p) 
U - sup - 

V {<p m 

For a qubit, the norm can be given by 



IA 



crio 



(13) 



(14) 



For a given system, the norm ||cr|L increase with time, 
reflecting the decoherence of the system. However, in 
general it is oscillated at the system's internal frequency. 
Thus, the decohering effect of the bath is better quanti- 
fied by the maximal operator norm, D (t) . This norm is 
also defined a measure for characterizing the short-time 
decoherence. 



D(T)=sup(||«7(r,p(0))y. 

P(0) 



(15) 



In the following, we will calculate the \\cr\\x then the de- 
coherence D of the JCQ. 



III. DECOHERENCE OF JCQ OPERATIONS 

As shown in above subsection, to calculate the deco- 
herence of the JCQ we must study its evolution under the 
interaction of the qubit with its environment. Suppose 
the initial state of the system be R (0) = p (0) ® 0, where 



2 



p (0) is the initial state of JCQ and 9 is the initial state of 
the environment, which is the product of the bath modes 
density matrices 9k- In the initial states, each bath mode 
k is assumed to be thermalized, namely, 

(16) 



where (3 = X/kT, k is the Boltzmann constant. So the 
evolution operator may be 

U - e -iHr/H _ e -i(H s +H I + H B )T/h_ ,yjs 

In the following we set t = t/H. Due to non-conservation 
of H s in this system, the evolution operator cannot be 
in a general way expressed as e - lH st e -i{Hi+H B )t_ g u |. m 
the sort-time approximation, the operator can be approx- 
imately expressed as [l4| 



U 



e -iH 3 t/2 e -i(H I+ H B )t e -iH 3 t/2 + o(t 3^ ^ 

It has been proved that the expression is accurate enough 
for the time being short to the characteristic time. So the 
density matrix elements of the reduced density matrix 
p (t) in the basis of operator H s can be expressed as 

e iH,t/2 e i(H I+ H B )t e iH,t/2 _ (1Q) 

By use of the completeness relation E |-) (-| = 1, we have 



±jH,t/2 



E e 

3=0,1 



where 



and 



Ao.i = ±- 



13 , 



±- 



Ej 



(20) 



(21) 



K) = _(| )-|1)), 



\<Pi) 



V2 



(|o) + |i». 



(22) 



Similarly, 

n=0,l 

(23) 

where Xo i = ±1 and \n) = |0) or |1) , are the eigenvalues 
and eigenstates of operator a z - So we have 

p mn (t) = Tr B (<P m \ E ^|^)(^| 

Q=0,1 

E e-^|^>(^| E |^><^| 

13=0,1 P, 9=0,1 

"pgr V"/ 1 1 1 " K z__/ " l^/i/ 



P M (0)I1^ E e' tA "|^>(^ 

k |U=0,1 



^ e *(xvE*J*+E*«*)* 

s=0,l 

E e^KX^IVn}, 
y=0,l 



(24) 



namely, 



a,/3,^,^,p,g,/x,z^— 0,1 

;>Iv„)<>„Iv„)tv^. ,-»(*E**+E»"0* 



Pp, (0)ll^e 



(25) 



In the Ea. f25p the Tr B term is same as that in Ref.Q- 
Enlightened by Privman's works we can easily obtain 

Pmn® = E e *(V^-A.-\0 

a./3,£,^,p,q,/i,f--0,l 

(Vml <Pa) ('Pal (£1 P/j) (^| <Pp) 
(<Pq\ Vp) \Vu| (?l Vv) (^1 <Pn) 

P P9 (0)e- B2 Wfe-^) 2 / 4 - iC? W(^-^)(26) 



Here, 



B *( t ) = sE^sin^-f coth^ . 

C(t) = X)^-(wfc*-sinw fc t), (27) 

and |0 , |?) € {|0) , |1» . S 2 (t) will affect the decoher- 
ence but C (i) will not it quantifies purely a shift of the 
energy levels of the qubit system. When the summation 
in Eq. I|27|l is converted to integration in the limit of infi- 
nite number of the bath modes, one has 



COht 



B 2 (t) = 8 / dioW (w) g (w) 2 u~ 2 sin 2 — coth — , (28) 



uit 

T 



(3oj 



for the real g (w) . By using of Eq. (|11|) we can yield a 
good qualitative estimate of the relaxation behavior 0, 
\LT\. From Ref.jl^] we know the condition of Eq.© is 
u) c h ^> Ej and uj c h ^> k B T. It is big enough for uj c to 
take uj c h — 200/iev. Thus, in the following numerical 
simulations we set the cutoff frequency uj c = 200 (the 
unit is the reciprocal of t's). Evaluating Ea. (|26|) we can 
obtain the evolution of the density matrix elements of 
the deduce density matrix p as 



Pw (t) 
Pu (*) 



2Pio( 
2 P °° ( 



e ~ B (*) + e itE J 



itEj-B 2 (t] 



-B 2 (t) 



->Pll 



1 



where p xx = p n (0), p 10 — p w (0) . The evolution of the 
closed qubit is p\ x (i) = p tl , and p\ (t) = p w e ltEj . So 
we have 



CT10 (t) 
Oil (t) 



Pw[l-e-» ^1(1 



-s 2 (*) 



Pll)- 



(30) 
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Then, we have 

{(Poo-Pii) 2 +4|p 10 | 2 sin 2 ^| 2 (31) 

The result is similar but not same to the norm in |7| 
where H s = — ^a z and A s = a x . 

In the following, we numerically analyze the decoher- 
ence. In the following calculation, three pure initial states 
are chosen, they are \<p) = (1,0) , (corresponding to 

point); \tp) 1 — (^p, , (corresponding to above line); 

1^)2 = (yT' V^) ' ( corres P oncun g to below line), [be- 
cause it has been shown that evaluation of the supremum 
over the initial density operators in order to find D, see 
Eq.(|15|l one can do over only pure-state density opera- 
tors l3]- We choose Ej = 51.8/xe« according to [l5| . 
and T — 30mK. In the model the typical impedance 
of the control line is R ~ 50S1. Since Rq ss G.bkil, 
Y. Maklin et al. have suggested a value n w 10~ 6 for 
numerical simulations. We plot the norms ||cr|K versus 
time t in Fig.l. 



Figl.a, Figl.b 



lv)oo (vl ' II^WIIa ' s ^ ne max i m u m and it equals to 
D = 5 fl — e-^WJ (plotted by points in the Figs.). We 

denote the low decoherence (D < 10~ 4 ) time t m . From 
Fig.la we obtain r ld w 7.5 x 1CT 2 x 6.582 x 10~ 10 s = 
49.4ps. We denote the elementary gate operation time, 
the characteristic time t 9 . In this case, t 9 = H/Ej w 
12.7ps. It is shown that the low decoherence time is larger 
than the single gate operation time, namely, r ld > t 9 . 
The fact means that within the whole time of elementary 
gate operation, D < 10~ 4 . Theoretically, the design can 
satisfy the DiVincenzo low decoherence criterion. A fur- 
ther study shows that when the temperature decreased 
and the Josephson energy Ej increased the qubit can be 
improved. It is also shown that the decoherence increase 
with the increasing of dimensionless strength of the dis- 
sipation rj. 

IV. CONCLUSIONS 

In this paper we investigated the short-time decoher- 
ence of the JCQ. We show schematically the behavior 
°f \\ a WIIa- It i s sri0wn that for a JCQ the decoherence 
derive from the dissipation is small enough according to 
the DiVincenzo criterion. It has been shown that the 
decoherence of JCQ derive not only from the dissipation 
but also from the quantum leakage. Fazio et al. pointed 
that the decoherence from quantum leakage is also not 
serious in JCQ model [l^. The two aspects information 
shows that the JCQ may be a good candidate of qubit 
for quantum computation. 



Fig.l: Norms versus time t, where the points 

and lines correspond to different initial states (see 
the text), Ej = 5l.8fj.ev, T = 30mk and ij = 10 -6 . 
The unit of the time in the Figs, is 6.582 x 10~ 10 s. 

It is shown that when the initial state is p (0) = 
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